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We evaluate the production cross section for direct Jj iff and T integrated in Pj for various collision 
pH energies in the QCD-based Colour-Singlet Model (CSM). We consider the LO contribution from 

gluon fusion whose Pj -integrated cross section shows a very good agreement with the Tevatron 
Q_i and LHC data, both for Jfifr and T. The rapidity distribution of this yield is evaluated in the central 

region relevant for the ATLAS and CMS detectors, as well as in the more forward region relevant 

for the ALICE and LHCb detectors. The results obtained here are compatible with those of other 
y—\ approaches within the range of the theoretical uncertainties which are admittedly very large. This 

1^ suggests that the "mere" measurements of the yield at the LHC will not help disentangle between 

t— \ the different possible quarkonium production mechanisms. Yet, the comparison with the first 

LHC results by ALICE, ATLAS, CMS and LHCb confirms that the CSM correctly accounts for 

the P T -integrated yield at ^[s = 7 TeV. 
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1. Introduction 

In 2007, the first evaluations of QCD corrections to quarkonium-production rates at hadron 
colliders became available. It is now widely accepted that a s and ar s corrections to the CSM [1] 
are fundamental for understanding the Pj spectrum of J /iff and T produced in high-energy hadron 
collisions [2, 3, 4, 5, 6, 7], while the difficulties of predicting these observables had been initially at- 
tributed to non-perturbative effects associated with channels in which the heavy quark and antiquark 
are produced in a colour-octet state [8, 9, 10, 11]. Further, the effect of QCD corrections is also 
manifest in the polarisation predictions. While the J/if/ and T produced inclusively or in association 
with a photon are predicted to be transversally polarised at LO, it has been recently emphasised that 
their polarisation at NLO is increasingly longitudinal when Pj gets larger [4, 5, 12, 13, 14]. 

In a recent work [15], we have also shown that hard subprocesses based on colour singlet QQ 
configurations alone are sufficient to account for the observed magnitude of the P^-integrated cross 
section. In particular, the predictions for the J/ip yield at LO [1] and NLO [2, 3, 4] accuracy are 
both compatible with the measurements by the PHENIX collaboration at RHIC [16] within the 
present uncertainties. This also pointed at a reduced impact of the ^-channel cut contributions [17] 
as well as of some specific colour-octet mediated channels relevant for the low Pj region Csi^ and 
3 pjy 8] ). The latter are anyway very strongly constrained by very important recent e + e~~ analyses [18] 
which leave in some cases no room at all for colour octets of any kind. 

The compatibility between the LO and NLO yields provided some indications that the com- 
putations are carried in a proper perturbative regime, at least at RHIC energies. The agreement 
with the data is improved when hard subprocesses involving the charm-quark distribution of the 
colliding protons are taken into consideration. These constitute part of the LO (a 3 s ) rate and can be 
responsible for a significant fraction of the observed yield [15, 20]. 

We proceed here to the evaluation the P^-integrated yield at higher energies both in the central 
and forward rapidity regions. At large energies, our study shows that the theoretical uncertainties 
on the J/ip yield for become very large -close to one decade- reminiscent of the case of total 
charm production [19]. Both for the J/ip and the T, we find a very good agreement with the CDF 
measurements [21] and the first LHC ones. Finally, we shortly discuss the impact of higher QCD 
corrections and the comparison with other approaches. 

2. Total J/ip and Y cross section at the LHC 1 

The Pj integrated cross sections obtained here have been evaluated along the same lines as 
our previous study [15]. The uncertainty bands have been evaluated following exactly the same 
procedure using the same values for m c (mb), /ur and //f. 

In Fig. 1, we show dcr C jJ^ ct / dy\ y= o x Br as function of -y^ from 200 GeV up to 14 TeV com- 
pared to the PHENIX [16], the CDF [21] and ALICE [22] data multiplied by the direct fraction 23 

'Note that we have not depicted the systematic uncertainties attached the unknown J /if/ and Y polarisation produced 
at the LHC . They can be as high as 50 % at very low Pj for extreme configurations. 

2 We employ this makeshift, although the^c yield can now be computed at NLO accuracy [23], since one cannot 
extend these computations to low Pj. 

3 Note that the measurement of the prompt yield by CDF went only down to Pj = 1.25 GeV. We have assumed a 
fraction of non-prompt Jf\fj of 10% below. We have assumed the same fraction at ^fs = 7 TeV for the ALICE data. 
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Figure 1: drf'™ ' / dy\ y=0 x Br from gg fusion in Figure 2: daf rect /dy\ y=0 x Br from gg fusion 

pp collisions for -\fs from 200 GeV up to 14 TeV in pp collisions for from 200 GeV up to 14 

compared to the PHENIX [16], CDF [21] and AL- TeV compared to the STAR [27], CDF [28] and 

ICE [22] data multiplied by the direct fraction (see CMS [29] data multiplied by the direct fraction (see 

text). text). 



measured by CDF at = L8 TeV [24]. We have found a good agreement. At larger energies, 
these results at 7 TeV (100 to 800 nb) and at 14 TeV (200 to 1400 nb) are in the same range as 
those of the Colour Evaporation Model [25] with central (upper) values of 140 nb (400 nb) at 7 
TeV and 200 nb (550 nb) at 14 TeV. They are also compatible with the results of the "gluon tower 
model" (GTM) [26], 300 nb at 7 TeV and 480 nb at 14 TeV, which takes into account some NNLO 
contributions shown to be enhanced by log(s). Quoting the authors [26], "the expected accuracy 
of the prediction is about a factor of 2-3 in either direction or even worse." In Fig. 2, we show the 
same predictions for direct Y for from 200 GeV up to 14 TeV compared to the STAR [27], 
CDF [28] and CMS [29] data multiplied by the direct fraction measured by CDF [30] at = 1.8 
TeV. While there seems to be some tension between the CSM predictions and the data at -y^ = 200 
GeV, the agreement is very good with the preliminary CMS data. 
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Figure 3: dcrjj™ c, /dy x Br from gg fusion LO con- 



Figure 4: d<7^ rect /dy x Br from LO CSM via 



tributions in pp collisions at yfs — 7 TeV com- fusion in pp collisions at V* = 7 TeV compared to 

pared to ALICE [22], ATLAS [32], CMS [31] and the preliminary CMS results [29] multiplied by the 

LHCb [33] results multiplied by the direct fraction direct fraction (see text), 
(and the prompt fraction (10 %) if applicable). 

In Fig. 3, we show the cross section differential in rapidity at = 7 TeV compared to the 
ALICE [22], ATLAS [32], CMS [31] and LHCb [33] results multiplied by the direct fraction (and 
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the prompt fraction (10 %) if applicable) as done for Fig. 1. The region y < is not plotted since 
it does not contain any additional physical information. Similarly, experimental measurements for 
y > or y < are physically strictly equivalent. In Fig. 4, we show the same plot for the direct T 
yield along the only public results available so far, from CMS [29], multiplied by the direct fraction 
from CDF. In both cases, the agreement is found to be very good. 

3. Discussion and conclusion 

Let us now briefly discuss the expectations for the results when QCD corrections are taken into 
account. First, we would like to stress that, although NLO results [2] are perfectly well behaved in 
nearly all of the phase-space region at RHIC energies [15], it does not seem to be so for larger s for 
the J /if/. One observes that the region where the differential cross section in Pj and/or y is negative 
(i.e. very low Pj and large y) widens for increasing s. Negative differential cross sections at low Pj 
are a known issue. Nonetheless, for yfs above a couple of TeV, and for some (common) choices of 
/up and hr, the Py-integrated "yield" happens to become negative, even in the central region. This 
can of course be explained by a larger contribution from the virtual corrections at -which can be 
negative- compared to the real emission contributions -which are positive-. Naturally, such results 
cannot be compared to experimental ones. This also points at virtual NNLO contributions at low 
Pj which are likely large but which are not presently known. Yet, as already mentioned, specific 
NNLO contributions were shown [26] to be enhanced by log(s). We also note that this issue of 
negative cross sections at LHC energies does not seem as critical for the T, whose production is 
initiated by gluons with larger xbj- 

As we have discussed above, one may try to compare the LO CSM with other theoretical 
approaches such as the CEM [25] and the GTM [26]. They all qualitatively agree, as well as with 
existing measurements. For all approaches, one expects a significant spread -up to a factor of ten - 
of the results when the scales and the mass are varied. 

Owing to these uncertainties, it will be difficult to discriminate between different mecha- 
nisms [8, 34] by only relying on the yield integrated in Pj and even, to a less extent, on its Pj 
dependent counterpart. This is a clear motivation to study at the LHC other observables related to 
the production of J /if/ such as its production in association with a single charm (or lepton) [15], 
with a prompt isolated photon [12, 13] or even with a pair of cc [3]. 

Yet, the good agreement at the level of the Pt -integrated cross section between the CSM 
and the available measurements from - 200 GeV to = 7 TeV confirms the little -if not 
negligible- impact of colour-octet contributions at low Pj, at least for the J /if/, in accordance to 
recent e + e~ analyses [18]. 
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